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Abstract 

It is shown that twist-3 fragmentation functions in the large-A'"c hmit of QCD obey the 
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^ . DGLAP evolution equations. The anomalous dimensions are found explicitly. The Gribov- 

t^^ ' Lipatov reciprocity is violated at the twist-3 level in the leading logarithmic approximation. 
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1. A renewed interest in the high energy spin physics in the last years has been concentrated 
on the transverse spin phenomena in hard processes. It revives many ideas developed over 
a decade ago. In particular, the notion of the twist-2 transversity distribution hi{x), first 
mentioned by Ralston and Soper |]I[ has been reinvented as well as its twist-3 counterparts 
have been addressed 0). Due to chirality conservation hi{x) cannot appear in the inclusive 
deep inelastic scattering (DIS) but it can be measured, for instance, in Drell-Yan reactions 
through the collision of the transversely polarized hadrons |jl|, |^ and in semi-inclusive pion 
production in DIS on the nucleon 0. In the last case, it enters into the cross section as 
a leading contribution together with the twist-3 chiral-odd spin-independent fragmentation 
function T{() |^. It is well known that there is considerable difference between the structure 
and fragmentation functions. Namely, the moments of the former are expressed in terms of 
reduced matrix elements of the tower of the local operators of definite twist. This property is 
established by exploiting the Wilson operator product expansion for inclusive DIS. Although 
the light-cone expansion for the fragmentation processes is similar to DIS, the moments of the 
corresponding functions are not related to any short-distance limit. As a substitute for the 
local operators come the Mueller's time-like cut vertices [Q which are essentially nonlocal in 
the coordinate space so that the analogy to the operator language is only useful mnemonic. 
Nevertheless, it is possible to give the definition of the twist for them in a broader sense without 
appealing to the concept of local operators. A thorough discussion of this issue can be found 
inRef. [|. 

The Q^-dependence of the twist-3 distributions has extensively been discussed in the litera- 
ture 1^, 1^, H, H, |T0[- As has been first observed in Ref. 0, in multicolour QCD {N^ -^ oo) the 
fiavour nonsinglet twist-3 part of the structure function (72 (x) obeys a simple DGLAP evolution 
equation, as in the case of the twist-2 operators, and the corresponding anomalous dimensions 
are known analytically. The same pattern is followed by the chiral-odd distributions [^ p!0|| . 



Since the twist-3 fragmentation functions enter into several cross sections on the same footing 
as the distributions, their scale dependence is of great interest. Apart from significance for 
phenomenology, it is important for theoretical reasons: while it is know that in the leading 
order of the coupling constant the splitting functions for the twist-2 fragmentation functions 
can be found from the corresponding space-like quantities via the Gribov-Lipatov reciprocity 



relation [jTT|, no such equality is known for higher twists. 

2. As distinguished from the leading twist evolution, the twist-3 two-quark fragmentation 
functions receive contribution from the quark-gluon correlators even in the limit of asympto- 
tically large momentum transfer. To solve the problem, one should correctly account for the 
mixing of correlators of the same twist and quantum numbers in the course of renormalization. 



Recently, we have addressed ourselves |T2[ to the question of the Q^-dependence of the nonpo- 
larized chiral-odd (NCO) fragmentation function 1{C) which is given in QCD by the Fourier 
transform along the null-plane of the appropriate matrix element of the parton field correlator 
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AO = \j'^e^'Hm{>^n)\KX){h,Xmm. (1) 

The summation over X is implicit and covers all possible hadronic final states populated by 
the quark fragmentation. We use the light-cone gauge B^ = 0, otherwise a link factor should 
be inserted in between the quark fields to maintain the gauge invariance. By exploiting the 
equation of motion for the Heisenberg fermion field operator we can express T{() in terms of 
the three-parton correlator Z{(', () which explicitly involves the gluon field, namely (neglecting 
the quark mass) 

i{0 = JdCz{C,0, (2) 

where we have introduced the C-even quantity 



z{cx) = ^[z^'\c'x) + [z^'\cx)r 



with 
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1 f d\dfi^^i^^_i^^, 



(0|g7+7p^^(An)|/i,X)(/i,X|^(0)R^(/xn)|0). 
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As we have observed in Ref. ||12[, in multicolour limit, i.e. neglecting the terms 0{1/N^), 
an additional three-parton correlator of the type {0\ipip\h, X) {h, X\B^\0) , which appears only 
through the radiative corrections, decouples from the evolution equation for Z{(', () and the 
latter becomes homogeneous provided we discard the quark-mass effects. Thus, in the large- iVc 
limit the RG equation takes the form 

and the evolution kernel is given by the following expression: 
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Here we have used the shorthand notation Gj^j^ for the following step function 

Inspired by our knowledge acquired from the study of the twist-3 structure functions [|, |, 
p!0[ , where the solution of the asymptotic equations was given by the convolution of the three- 
particle correlator with a certain weight function that is essentially the same as entering into 
the equation that gives rise to the dynamic twist-3 contribution to the two-quark distribution 
at the tree level, we are able to check that Eq. (0) satisfies the ladder-type evolution equation 
with the following splitting function: 



^JdCic{z,z'xx') 



i 1-i 



^^-iH^-l 



Thus, for the moments we obtain the following solution of the RG equation {Q > Q, 
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and the corresponding anomalous dimensions equal 

^^°7n = N, {-2^(r^ " 1) - 27i. " ^^ + ^} ' (1°) 

as usual /3o = |A^/ — ^Ca- 

As we have previously mentioned, there exists an equation which states that in the leading 
log approximation the time-like (TL) and space-like (SL) kernels corresponding to the twist-2 
parton densities are directly related, in the physical regions of the corresponding channels, by 
the Gribov-Lipatov equation P^^{x) = P"'"'"(l/x), or in terms of the corresponding anomalous 
dimensions it looks like 7J+2 — 7n^- Comparing the result given by Eq. (|T0|) with the large- 
Nc anomalous dimensions known in the literature for the chiral-odd distribution e{x) |^, |TU| 



we see that there is no universality of the corresponding twist-3 evolution kernels, i.e. the 
Gribov-Lipatov reciprocity is violated. 

3. Now we can proceed further and demonstrate that the evolution kernels for the time- 
like two-quark densities can directly be found from their space-like analogues by exploiting the 
particular form of the evolution equations given in Ref. [|10|. Since the analytic structure of 



the uncut diagram (see fig. |l|) is completely characterized by the integral representation of the 
G-function given by Eq. (^), we can just take its particular discontinuities, using the usual 
Cutkosky rules supplied with appropriate theta-function specifying the positivity of the energy 
flow from the right- to the left-hand side of the cut, in order to obtain the corresponding time- 
like kernel. Since the observed particle is always in the final state for the fragmentation process. 




Figure 1: One-loop ladder-type diagram for the two-particle evolution kernel. 

we are restricted to the single cutQ across the horizontal rank of the ladder diagram. Namely, 
using the integral representation of the corresponding step function (0), we have 

^0/ _B\=r— i disc e{x-(3) 

^^^"^'^ ^' J-oo2m[ax-l+t0][a{x-f3)~l + t0]^ (3 ' ^ ' 

The self-energy insertions are not affected by the cut since it does not cross the corresponding 
lines. Taking into account different kinematic definitions of the correlation functions in the 
space- and time-like regions [|T^, we are able to find the kernels. It is easy to verify that the 
evolution kernels constructed for the time-like twist-2 cut vertices using this recipe coincide 
with the known results. In the same way, we may obtain the above equation (H) from Eq. (90) 
ofRef. |ig. 

Since there exist fragmentation functions corresponding to each distribution, apart from 
the specific ones appearing from the final state interaction, we are in a position to find large- 
Nc anomalous dimensions, which govern their Q^-dependence, from the results of Refs. [|, 
^ 1^, Q. Namely, the genuine twist-3 contributions 7i^~^ and Q^~^ to the corresponding 



fragmentation functions 

ndO = \ /^e^^^(0Ka+_75^(An)|/.,X)(/i,X|^(0)|0), (12) 

4 J /vr 

Gt{0 = 7 f^e^'HOb±l,H>^n)\h,X){h,X\m\0)- (13) 

4 J ZTT 

after subtracting out the twist-2 piece p obey the evolution equation (H) with the following 
anomalous dimensions: 

''''''in = iV,|-2^(n-l)-27i, + ^ + ^}, (14) 

''^''7n = N,[-2^Pin-l)-2^E--^ + ^Y (15) 

Of course, it is a trivial task to invert the moments and to find the DGLAP kernels themselves. 



^For a given graph the possible cuts correspond to the possible final states. 



To summarize, we have found that in the multicolour limit of QCD the twist-3 fragmentation 
functions obey the ladder-type evolution equations and corresponding anomalous dimensions 
are known analytically. This gives a possibility to analyze experimental data when they become 
available. The Gribov-Lipatov reciprocity is not the property of the twist-3 distributions but 
rather it is strongly violated already in the leading logarithmic approximation. 
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